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The radiative instability of the relativistic electron beam in a periodic dielectric- 

\ filled cylindrical waveguide is considered. The dependence of the beam instability 

. ■ increment on the radiated wave frequency near the region of dispersion equation 

roots degeneracy is studied. It is shown that sharp change in the instability of the 

beam under the conditions of two-wave diffraction in Compton generation regime, 

making the increment proportional to the fourth root from the beam density p^/^ 

in contrast to conventional law p^^'^, brings radiation generation (amplification) in 

Q \ the considered system to be essentially improved in comparison with conventional 

O 

I devices (BWO, TWT, FEL etc). Numerical calculations of the instability increment 

c/j ' 

O ■ for various parameters of the system are performed. 

I. INTRODUCTION 

> 

cn 

■ Currently, there is a large number of generators and amplifiers of electromagnetic radi- 

oo : 

CN ■ ation (from microwave to optical wavelengths range) based on electron beams, for example 

^ \ traveling- wave tubes (TWT), backward wave oscillators (BWO), free electron lasers (FEL), 
ubitrons etc . It is known that any radiating system is characterized by its dispersion 
equation describing in the case of small perturbations the possible types of waves in the 
system. A detailed analysis of the properties of this dispersion equation jsl shows that the 
gain in the Compton regime (increment of electron beam instability) of the most commonly 
used generators (TWT's, BWO's, FEL's) isproportional to p^, where p is the density of 
the electron beam. However, in the papers 0, it is found that for electron beam moving 
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in spatially periodic medium under the conditions providing the coincidence of roots of the 
dispersion equation there is a new physical law: the increment of instability is proportional 
to p^/^'^+^\ where s is the number of additional waves appearing due to diffraction in the 
crystal. The analysis shows that this new law leads to reducing of the electron beam current 
density, needed to reach lasing threshold. This enables development of a new type of Free 
Electron Lasers called the Volume Free Electron Lasers (VFEL) . Due to a significant 
change in the threshold conditions, VFEL can provide a more efficient radiation process 
than conventional generators. 

Until recently the theoretical studies of the problem of beam instability in photonic 
crystals were carried out for the case of infinite in the transverse direction crystals. However, 
in many cases the mode structure of the electromagnetic field in laser cavities can not be 
neglected (for example, when the lasing is performed in the microwave range). For the first 
time the process of VFEL lasing in photonic crystal finite in transverse direction is considered 
in 8] . In js I it is also pointed out that the four-fold degeneracy of the dispersion equation 
roots (when increment is ~ p^^^) is possible in the one-dimensional photonic crystal when 
the condition Eq > 1 (in the conventional generators - BWO, TWT etc. Eq is usually equal 
to 1) is fullfilled. In this paper we consider the simplest example of such one-dimensional 
photonic crystal - a cylindrical waveguide with a periodic dielectric filling. 



II. BASIC FORMULAS 



The system of equations describing the interaction of an electron beam with an electro- 
magnetic wave in a waveguide can be obtained from the Maxwell and electron movement 
equations. We rewrite the Maxwell equations as follows 

rot rot ii/(r, a; ) -E\r ,uj)E\r ,uj) = — —j[r,u), (1) 

div E{f,u)E{f,u) = 4:71 p{f,uj), (2) 

iujpif.uj) — div j(r,a;) = 0, (3) 

where E{f,uj) = J E{r^t)e^^^dt is the Fourier transformation of the electric field E{'r,t); 
E{f,uj) is the dielectric permittivity of the medium filling waveguide; j{f,u) and p{f,u) are 
Fourier transformations of the current density j(r, t) and electric charge density of the beam 
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p{r,t), respectively. The quantities j{r,t) and p{r,t) can be expressed in following way 

j(f, t) = eY, Va{t)S{r~ f^{t)), (4) 

a 

p{f,t) = eJ2^ir-f^{t)), (5) 

a 

where fa{t), Vait) are ath electron radius- vector and velocity; the sum in dl])-© is over all 

electrons in the beam. 

The electron movement equations can be written in the form 

dvait) e 
dt m'-f 

where E{fa{t),t) and H{fa{t),t) are the electric and magnetic field of the electromagnetic 
wave in the point ra{t) at the time moment t, 7 = (1 — 

Let the z coordinate axis coincide with waveguide axis. We also suppose that electron 
beam is "cold" (velocity spread of the electrons can be neglected) and initial electrons ve- 
locities are directed along 2;-axis {u = uCz] {gx, Cy, e^) are unit vectors of coordinate system). 
The dielectric permittivity inside the waveguide is e{f,u) = EqIp) + x(r, w), where x is 
periodic function of z: x(r, w) = Yl Xt{p)g''^^', Xo = ^0 — 1- With the help of and 
supposing that |x(r, a;)| ^ 1 we can rewrite ([1]) in the following way 

- V^Eif, u) - — £(r, u)E{f, u)^ — [ j(r, u) + ^V(Vj(r, u)) . (7) 



The field E{f,u) can be decomposed in terms of the waveguide eigenfunctions 

E{f, uj) = ^Y. f ^n{h)^n{p, h)e'''''dh, (8) 

n 

—* 

where an{kz) are the expansion coeffitients (arnplitudes), Yn{p, kz) and k„ are the waveguide 
eigenfunctions and corresponding eigenvalues |9l-lll|. 

The beam current appearing on the right-hand side of ([7]) is a complicated function of 
the field E. To study the problem of the system instability, it is sufficient to consider the 
system in the approximation linear over perturbation, i.e., one can expand the expressions 
for j(r, cj) over the field amplitude E{r,u): j = jo + 6j, where jo is the beam current not 
perturbated by the radiated field, 6j ~ E{f, u) is the beam current induced by the radiated 
field. With the help of (j4]) one can find 

6j{k, w) = e e"*''''"" < 6va{oj - ku) H (k5va{oo - ku)) \ . (9) 

I u -ku ^ ^ ) 
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The quantity Sva{u}) can be obtained from movement equations: 

= / AT^e { —E{k\ u + k'u) + 

mu-f J {27ry \uj + k'u \uj + k'u 



(10) 

After substituting ffTOjl into ([9]) in the expression for the current density appear the sum 
^gi(fc'-fc)r,o_ Let us average this sum over distribution of the particles in the beam 

a 

^gi(fc'-fc)r,o ^ ^(j:^ _ . {2Txfno5{k, - k'^), (11) 

a 

where ^{k± — k'j_) = Jge~^^''^~'^'^^^(p{p)(Pp, ^ Jg(p{p}cPp = 1, 5 is cross-section of 
the waveguide, Uq is the electron density of the beam, k^ is the longitudinal component of 
the wave vector, the function (p{p) describes the distribution of the electrons in the beam 
cross-section. 

The expressions (17|)- (|TT]) allow us to write the following system of equations for amplitudes 
0'm{kz) (this system is similar to that describing the multiwave dynamical diffraction in 
crystals): 

( h ^ ] ^mnO"n{kz), (12) 



where uf = "° is the Langmuir frequency of the beam, effective susceptibility xTpr is 



XTPfik,, k,-T)= j Y:,{p, k,)XrYn{p, k, - T)(fp. (13) 



Coeffitients Amr, in flT2l) are 



X k,){u- -^fc) (-ic'ix V'j. + [aVrAf', t)) . (14) 
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III. INCREMENT OF ELECTRON BEAM INSTABILITY 



Let us consider the generation on the most commonly used -Eoi-niode of circular waveg- 



uide. The vector eigenfunction Yi{p, kz) corresponding to -Eoi-mode can be written as [11| 



^1 = T / — V { /,o ^ Ji{p'^)^r + Joipi^)ez \ , (15) 

where i^oi ~ 2.4048, n = uqi/R, R is the waveguide radius, Jo and Ji are Bessel functions. 
The effective susceptibility according to (fT3|) is 

Xe//(A:., r) = Xrj Y*{p, kz)Up, k^ - T)d'p = Xr (^1 - ^2^) • (16) 
The equations (fT2|l for the i?oi-mode can be written as 



[CO 



kzuf {kl - k%)a{kz) -klY\ Xeffikz, T)a{kz - r) = -^Aiia(fc^), (17) 



where /c^o = y — and /cq = u)/c. We assume that unperturbed beam occupies the 
entire cross section of the waveguide and (p{p) = 1. Then the calculation of the An in 
accordance with f[T^ gives 

Let us consider firstly the case of waveguide with homogeneous filling {xr = 0). Ap- 
parently in this case the diffracted wave is absent and system (1171) is reduced to a single 
equation. Suppose that right-hand side of f[T7|) is small. Since the nonlinearity is insignifi- 
cant, let us consider as the zero approximation the spectrum of the waves of equation f fT7|l 
with zero right-hand side. So, the solution of (|T7|) is kzi = k'^^ + Skzi, where 6kzi is to be 



found, and fc^i is the solution of f|T7|) with zero right-hand side k'^^^ = kzo = \ ^£^o — 



Repeating exactly the arguments, for example js], we arrive the following expression for 
increment 



,2 



Now let us consider the case Xr 7^ 0, Xt ^ 1- Apparently when the diffraction conditions 
are not fullfilled we get the same expression 0191) for increment. We therefore assume that 
the conditions of two-wave dynamical diffraction in a waveguide are realized, i.e. the wave 
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amplitude a{kz + r) is comparable with the amplitude a{kz). In this case the dispersion 
equation has the next form: 

- k^uf {{kl - kl,){kl - kl) - klxefAk^r. r)Xeff{k., -t)} = "^^lllC " k'J ■ (20) 
The solution of this equation is kz2 = k'^^ + ^kz2, where 



4'''^ =ho\l- \aB(3 ± ^iaB(3)^ + ^^(3 } , (21) 



= Xeffikzo + T, T)xeff{kzo, -t) , k^Q = J - K^, 73 = /? = is the diffraction 



asymmetry factor, Ob = ''^^f'a^'^'*'^ is the off-Bragg parameter (as = when the Bragg 



zO 



condition of diffraction is exactly fulfilled), t = 27t/D, D is the waveguide period. The 
calculation gives the following expression for increment Imkz in the vicinity of the roots 
degeneracy point 

Im kz2 = Im 6kz2 ^ (^^i^^n) ~ P'^'- (22) 
From ([H]) and ([22]) we have 

« » 1 (23) 

because uf <^ w^. An <^ /cq. This means under diffraction conditions the gain increases. As 
a result, the threshold current density in two- wave diffraction case is lower (j ~ (fcL)«(fcx l)^" ' 
where L is the interaction length). 



IV. NUMERICAL SOLUTIONS 



To be able to experimentally observe this effect it is necessary to know how accurately 
the diffraction condition must be satisfied. To study this, the equations fl20|) were solved 
numerically. Now we demonstrate the results of this calculation by the example of one of 
the systems studied (i? = 6 cm, D = 3.6 cm, eq = 1.23, Xt = 0.05). 

On the fig. [1] one can see the dependence of the increment of instability on the frequency 
near region of dispersion roots degeneracy (4.06 GHz in our case). Let us explain how the 
curves in fig. [1] are obtained. Initially, we found the solutions of dispersion equation without 
beam ( 12T|) . They are shown on the right-hand side of the graph. Next, for each frequency 
fi from the selected range the set of equations fl2Ul) was solved and the imaginary part of 
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FIG. 1: The dependence of instability increment of electron beam on the frequency near the point 
of the roots degeneracy of dispersion equation. Beam current is 0.1 kA. Solid and dashed curves on 
the right-hand side of the graph are dispersion characteristics of the forward and backward waves 
in the waveguide, respectively; on the left-hand side are plotted corresponding increments. On the 
graph are also shown the beam lines (w = k^u) for the energies 600 keV and 1 MeV and the light 
line {io = kzc). 

the solution found Im kz was plotted. In the calculations we assumed that the beam current 
I is given; the beam energy E was chosen so that the Cherenkov synchronism occurred at 
a frequency /j. Calculations were performed for the two roots of the dispersion equation.^ 
For example, for a backward wave synchronism with the beam with the energy 600 keV 
occurs at a frequency of ~ 3.9 GHz; if we choose the beam current equal to 0.1 kA, the 
corresponding increment will be Imkz ~ O.OOGtt/D ^ 0.007 cm~^, which is shown in Fig. [H 

^ Far from the roots degeneracy point they correspond to the so-called forward and backward waves in the 
waveguide. But near the point of degeneracy such division, generally speaking, is incorrect because in 
this case "forward" and "backward" waves are coupled. Nevertheless, in the article we use the established 
terms "forward" and "backward" waves to refer to the roots of the dispersion equation, since we investigate 
the behavior of the system over a wide frequency range (far from the point of degeneration and in the 
immediate vicinity from it). 
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FIG. 2: The growth of the increment with beam current increase. On the graph are shown the 
increments for cases when the beam current is 0.1 A, 1 A, 10 A, 0.1 kA and 1 kA. 

A similar example is given on the fig. [T] for the forward wave. 

From Fig. [1] we can trace back the typical dependence of the increment value on the 
frequency near the point of the roots degeneracy for the systems studied. It is easy to see 
that the synchronism condition for the forward wave is possible only for frequencies above 
/o (in our case /o ~ 3.87 GHz), because velocity of the beam can not exceed the speed 
of light in vacuum c. Therefore, for the forward wave the increment begins to rise from 
zero at the frequency /o to a maximum which can be observed both at the point of the 
roots degeneracy, and at some distance from it (see fig. [2]). It is seen that in the vicinity 
of the degeneracy point the increment increases with the frequency faster than away from 
the diffraction conditions. For the backward wave a similar but much more pronounced 
picture can be seen: relatively slow increment growth away from the diffraction conditions 
is replaced by the rapid growth near the degeneracy point. In addition, for backward- wave 
the increment has the highest value exactly in the degeneracy point. 

Fig. |2] shows the increment of instability for different beam currents. For the forward 
wave the increase of the beam current shifts down the frequency at which the increment 
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reachs maximum. When the current rises, the increment maximum becomes wider and 
lower. Nevertheless, for the chosen parameters of the waveguide the gain for a backward 
wave in the vicinity of the degeneracy point can be several times (from 2 to 5) greater 
than that for the usual no-diffraction case. But to use this in practice the synchronism 
conditions must be satisfied very precisely. For example, when the beam current is 1 kA as 
shown in fig. [2] the imaginary part Im is reduced by half when the displacement from the 
point of degeneracy (4.06 GHz) down in frequency is only ~ 116 MHz (the peak FWHM 
is ~ 116 MHz). In practice, the control of oscillation frequency with such precision is not 
always possible due to many factors, such as the presence of an electron velocity spread 
in the beam. Nevertheless, one can find conditions (such as geometric dimensions of the 
system, values of Eq, Xt, the beam energy, etc.) under which it will be possible to perform 
the lasing near degeneracy point with sufficient accuracy. In particular, by use of fig. [T] it is 
easy to find that for selected geometry the spread in the electron velocities does not play a 
significant role due to high electron beam energy (for the detection of the effect is enough 
to have 5E < 30% at the energy E ~ 1.6 MeV). 



V. CONCLUSION 



We conclude with a remark on the above mentioned model of the one-dimensional crys- 
tal. This model is is somewhat idealized since in real systems it is impossible to use the 
waveguides with a solid dielectric filling. As a rule, in practice are applied waveguides with 



inner dielectric liner 



12| (partially filled waveguides) . However in considering of the partially 



filled waveguide all the above arguments remain valid. In fact, you only need to substitute 
in the above model other eigenfunctions and derive the expressions for coefficients Xeff and 
Ajnn- But the behavior of the instability increment near the degeneracy points remains the 
same. Moreover to estimate the increment value you can use the model of a waveguide with 
a solid filling by introducing the average over the cross section of the waveguide dielectric 
permittivity e{uj) and taking into account the fact that with a wave interacts effectively only 
part of the beam (located at distance of < A/37/(47r) from the surface of the dielectric). 

It is shown that for a one-dimensional photonic crystal the gain near region of dispersion 
equation roots degeneracy rapidly (several times) increases. With a reasonable choice of the 
waveguide (crystal) and electron beam parameters (in particular, at sufficiently high beam 
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energy) this effect is experimentally observable and can be used to improve the characteristics 
of generators and amplifiers of microwave radiation. 
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